Two-parameter Asymptotics 
in Magnetic Weyl Calculus 



Max Lein 



May 31, 2010 

Technische Universitat Munchen, Zentrum Mathematik, Boltzmannstrafse 3, 85747 Garching, Germany 

e-mail: lein@ma.tum.de 



Abstract 

This paper is concerned with small parameter asymptotics of magnetic quantum systems. In addition 
to a semiclassical parameter e, the case of small coupling A to the magnetic vector potential naturally 
occurs in this context. Magnetic Weyl calculus is adapted to incorporate both parameters, at least one 
of which needs to be small. Of particular interest is the expansion of the Weyl product which can 
be used to expand the product of operators in a small parameter, a technique which is prominent to 
obtain perturbation expansions. Three asymptotic expansions for the magnetic Weyl product of two 
Hormander class symbols are proven: (i) e < 1 and A< 1, (ii) e <C 1 and A = 1 as well as (hi) e = 1 
and A -C 1. Expansions (i) and (iii) are impossible to obtain with ordinary Weyl calculus. Furthermore, 
I relate results derived by ordinary Weyl calculus with those obtained with magnetic Weyl calculus by 
one- and two-parameter expansions. To show the power and versatility of magnetic Weyl calculus, I 
derive the semirelativistic Pauli equation as a scaling limit from the Dirac equation up to errors of 4th 
order in Vc 
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1 Introduction 

Quantum mechanical systems often contain small parameters that allow us to order terms by magnitude 
and importance. One prominent example are adiabatic systems where the fast degrees adjust 'instanta- 
neously' to the configuration of the slow degrees of freedom. Here, the small parameter quantifies the 
separation of slow and fast scales. Under certain conditions, effective dynamics may be derived which 
contain corrections order-by-order in the small parameter and one can bound the error. This effective 
hamiltonian may be the starting point for a semiclassical analysis: Egorov-type theorems compare the 
quantization of suitable classically evolved observables with the corresponding time-evolved quantum 
observables. 

A quantization procedure is a systematic way to associate operators to functions on symplectic mani- 
folds that has certain natural properties (e. g. linearity and compatibility with the involution, see Chap- 
ter 5 [|Wal08H for an overview). Mathematically speaking, we are interested in a functional calculus 
for non-commuting observables called position x and momentum E, on phase space T*R d = R d x R d 
endowed with the magnetic symplectic form. This is by no means the only interesting case, other ex- 
amples are spin systems HVGB89I IVGBS901I or quantization procedures on generic Poisson manifolds 
IIKon03l IWal08ll . Before we explain magnetic quantization in detail, we will briefly recall the non- 
magnetic case. 

Usual Weyl quantization Op E : / — > Op E (/) maps suitable functions / on phase space T*R d onto linear 
operators acting on (subspaces of) L 2 (_R d ) (see ||H6r79llFol89ll , for example). The index e indicates that 
the commutator of momentum Op E (£) and position Op E (x) is of order e, i [Op E (£ ( ), Op E (x ; )] = e5 i; . 
With this quantization procedure in hand, it turns out we can define a non-commutative product | e on 
phase space which emulates the operator product, 



If e <C 1, we can expand the Weyl product asymptotically in powers of e up to an arbitrarily small error 
which allows us to rewrite the operator product as an asymptotic series in e as well ||Fol89ll . This idea 
has been used to derive effective operators whose dynamics approximates the full perturbed dynamics, 
see, for instance, ||LW93| IPST03bl IPST03al ITeu03ll . Hence, from a computational point of view, an 
asymptotic expansion in a small parameter is a very desirable thing to have. 

Very often, the magnetic field B is the perturbation of the hamiltonian and usual (i. e. non-magnetic) 
Weyl calculus is not well-adapted to this situation. A magnetic Weyl calculus is needed and this paper 
is concerned with the derivation of three asymptotic expansions of the corresponding magnetic Weyl 
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Op E (/)Op E (g) = Op E (/tf E g) x £V Op E ((/tf E g) (n) ) +Op E (0(O). 
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product, e. g. with respect to a semiclassical parameter s, 

00 

n=0 

Let us introduce some notation first: assume we apply a magnetic field, thought of as a two-form 
B = (Bij)i<i;<d> whose components are smooth, bounded and have bounded derivatives to all orders, 
i. e. B (; e ^ < ^°°(IR d ), 1 < l,j < d. Then we will consider the quantization which takes the position and 
momentum vectors x and E, into the operators 

Q:=x (1.1) 
P A A := -ieV x - AA(Q) 

A is a vector potential which represents B, i. e. B; ; = d x Aj — d x .A h 1 < I, j < d, whose components will 
always be chosen to be smooth and polynomially bounded, A { e < ^ p c ^ 1 1 (M d ), 1 < Z < d. The first parameter 
e formally takes the role of ft and quantifies the non-commutativity of position and momentum. There 
are many systems (e. g. Born-Oppenheimer-type systems) where e is, physically speaking, not ft, but 
some other parameter that describes a separation of scales. The second parameter, A, is often physically 
equal to e /c where e is the charge quantum and c the speed of light. If A is taken to be small, this could 
either mean that we are interested in some non-relativistic limit, c >■ 1, or in the limit of small charge 
e<l. To be precise, small parameters must not have any units, so instead of A = l /c we should really 
use A = v o/ c where v is some characteristic velocity. 
The commutators of our basic building blocks are given by 

i[Qi,Q } ] = 0, iLPf^Qi] = e5 Xj , ^K,xi'Kx^ = -**By(Q)- 

We note that the commutator of the kinetic momentum operator depends on the magnetic field B and 
not on the specific choice of vector potential. If we translate these commutation relations to the classical 
framework where x corresponds to Q and E, to Pf x , we have to use the magnetic symplectic form 

d d 

oj xb = ^d£, A dx, - ±A V B l] {x)dx l A dx ; - 
(=i i,j=i 

on phase space T*R d which induces the magnetic Poisson bracket, 

d d 

{f>g}, B = E( V \S - V %g) - A 2 BijW %/ d ljg . 
i=i i,j=i 

Classically, this agrees with the recipe of minimal substitution: if we replace £ with £ — AA(x) and use 
the standard Poisson bracket {•, •}, we recover the magnetic Poisson bracket, 

{/(*, ? - AAM), g{x, I - AA(x))} = {/, g} XB {x, I - AA(x)). 

Quantum mechanically, these two points of view are no longer equivalent. Based on the magnetic sym- 
plectic form, Miiller was the first to define covariant magnetic Weyl calculus in a non-rigorous fashion 
||Mul99ll (the author thanks R. Littlejohn for this reference), although Luttinger has used it in proto- 
typical form as early as 1951 IILut51ll . The present paper relies upon earlier contributions which have 
put these ideas on a solid mathematical foundation HIMP071 IMP041 IMP051 IKO01I IKO041 IKO05II . Some 
notable results include a Calderon-Vaillancourt-type theorem (L 2 continuity of 5^° symbols), selfadjoint- 
ness of elliptic symbols on magnetic Sobolev spaces HIMP07II and a Beals-type criterion IIIMP08II . One 
missing ingredient is an asymptotic expansion of the magnetic product with respect to a small parameter. 



3 



Before we continue, we would like to elaborate on possible choices of scalings. If we rescale space by 
e via := e d/ ' 2 tp{ex), y> e L 2 (M d ), we can transform the observables (II. ID into 



Q E := ex 



(1.2) 



if 



-iV x - AA(Q e ). 



Mathematically, both scales are unitarily equivalent (see Appendix©. The decision which scale is 
deemed preferable is based on the physics of the problem. If we would like to emphasize the slow 
variation of the magnetic field (compared to other potentials), then the second choice is more natu- 
ral. The single-particle Schrodinger equation with periodic potential V r subjected to a slowly-varying 
electromagnetic field, a system which is described by the hamiltonian 

H = H~iV x - A(ex)) 2 + V r (x) + <*>(ex), 

falls in this category. We emphasize that all of our results hold in either scaling. In particular, the 
asymptotic expansion of the product is the same, independent of the choice of scaling (see Appendix |A| 
for details). As this paper was initially motivated by the problem above, we will use the adiabatic scaling 
given by equation HI. 2D . 

The fundamental building block of magnetic pseudodifferential calculus is the magnetic Weyl system 
that implements the commutation relations, 

W* x {X):=e- ia( - x ^KJ\ 

Its gauge-covariance leads to the gauge-covariance of magnetic pseudodifferential operators. Here 
(x,^) = X e T*R d is a point in phase space and cr{X,Y) := £ • y - x • rj, Y = (y,rj) e T*R d , is the 
(non-magnetic) symplectic form. In HMP04II it has been shown that this is a well-defined operator which 
acts on any <p e L 2 (M d ) by 



The appearance of the magnetic circulation V A along the line segment which connects x and y 



Ar A ([x,y]):=- ' 

E e 



A 



(1.3) 



stems from the use of magnetic translations. All proofs in [IMP07, MP041 IMP05II carry over to the present 
case via a simple scaling argument. The pseudodifferential operator associated to a Schwartz function 
/ e 5"(T*IR d ) is defined in terms of the symplectic Fourier transform & a f = &~ l f and the Weyl system: 



(27l) d 



dX(^fXX)Wf x {X):= — 



c\Xe iaixJi) f{X)Wf x (X) (1.4) 



All parameters are contained in the Weyl system Wf x {X~) = e 1<T(x ' (<J " n a"; if we had chosen the usual 
scaling, the formula would be the same, but Q E and n A A would have to be replaced by Q and Pf x from 
equation HI. ID . This definition can be extended to observables of Hormander symbol class m with weight 
p e [0,1] [M P04lliMP071] among others: 

5" p m := {/ G ^°°(T*]R d ) I Va,aeN d 3C aa > : \d x a d«f(x,Z)\ < C aa (£)"H a |p} 
The quantization of Hormander-class symbols / act on Schwartz functions <p e 5^(R d ) as 



(op^cn^xx): 



(27l) d 



drj e 



-i{y-xyT] p -iXY A e {ix,yYj 
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where the integral is interpreted as an oscillatory integral. If in addition the symbol / is elliptic of 
order m, for instance, this definition extends from xp e 5^(M d ) to functions in the (magnetic) Sobolev 
space H™(R d ). These results are well-known for standard Weyl calculus, see e. g. ||H6r72| IDH731 |H6r79l 
IH6r85l|Ste93l| . If we choose an equivalent gauge, i. e. A'(x) = A(x) + eV x xM for some % e < ^ p ~(K d ), 
then the magnetic Weyl quantization of / with respect to A' is related to that with respect to A by 
conjugating with e +ax ^°\ that is, magnetic quantization is covariant: 

Opf x eV ' x if) = e +a * CQ - 3 Op^(/)e- lA ^ Q ' ) 

Unless f is a polynomial of degree < 2 in momentum, regular Weyl quantization of minimally substituted 
symbols is not gauge covariant and does not coincide with magnetic quantization. This includes physically 
relevant examples such as \J m 2 + E 2 or band energy functions in solid state physics. 

The second major component in magnetic Weyl calculus is a product jj^,. Its form is shaped by the 
commutation relations of the fundamental observables as expressed by the composition law of the Weyl 
system, 

W A A (T) W A A (Z) = e'S CT(nz) e ->7r*(<Q ( .,Q £ + E y,Q ( .+ £ y+ ez >) W f x {Y + Z). (1.5) 

The magnetic contribution e - | 7 rB ( (Q " Q > +E:> '- Q » +e;> ' +E2:> ) is the exponential of the magnetic flux through the 
triangle with corners Q E ,Q £ + ey and Q e + ey + ez (see equation (I2.2D ) and as such depends only on the 
magnetic field and not on the choice of gauge. A scaled magnetic flux through triangles with different 
corners (equation (12.3D ) enters the integral formula for the magnetic Weyl product, 



i r 



(2ti) 



2d 



c\Y 



dz e+i(7 (x,Y+z) e q apw-ttrfcxovo f>;7)(Y) (J^gKz). (1.6) 



This expression can be derived from the composition law of the Weyl system and the magnetic Wigner 
transform (see Theorem 12. 10D . If e •< 1 and y and z have lengths of 6(X), then j B E {x,y,z) is of order 
6(e) and can be expanded in powers of e. The different asymptotic expansions with respect to e and A 
simultaneously, e or A are obtained by expanding exponential of the 'twister' 



e 



iT eJ Xx,Y,Z) ._ i^cj{Y,Z)-iXy B e {x,y,z)) 



in the appropriate way. One can easily distinguish purely magnetic and purely non-magnetic contribu- 
tions in the product. 

In the 'usual scaling' (equation HI. Il l), the Weyl system would essentially obey the same composition 
law and lead to the exact same expansion of the magnetic product (see Theorem lA.3l for details). 



1.1 Comparison with usual Weyl calculus 

Regular Weyl quantization has seen many applications to magnetic systems over the years, so one obvious 
question is how results would differ if magnetic Weyl calculus had been used instead. In ordinary Weyl 
calculus, a magnetic field is included by quantizing the minimally substituted symbol h o # M instead of 
h where -& M (x, := (jc, E, — AA(x)) . If h e is a Hormander symbol, then one typically assumes that 
the magnetic field B is such that there exists a vector potential whose components of A are ^ < ^ 00 (R d ) 
functions or that all derivatives are bounded, |3 x a A(x)| < C a for some C a > and all a e Nq, \a\ > 1. 
The latter is typically used to cover the case of constant magnetic field. Under these assumptions, 
Theorem |2 . 1 9 1 ensures the difference 

Op E (ho-& M )-Opf A (h) = ^(e 2 ) (1.7) 

is a pseudodifferntial operator which is the quantization of some g e 5^ m ~ 3 . The difference can be 
expressed as a power expansion in e. 
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The second important ingredient is a comparison of the two products: standard Weyl calculus suggests 
to multiply / o # M and g o # M where / e 5^ mi and g e 5^" 2 . Combining 

with equation Q1.7D , we also get that the two products differ by 0{e 2 ), 

Op E (/ °# M tU = Op^C/fl^g) + 0{e 2 ). 

Hence, first-order corrections (e. g. the subprincipal symbol of an effective hamiltonian) derived with 
usual Weyl calculus coincide with those obtained by magnetic Weyl calculus. However, the advantages 
of magnetic Weyl calculus are significant, most importantly do results extend to much more general 
magnetic fields, one only needs to assume that the components of B are in ^ < ^ 00 (ffi d ). The gauge- 
covariance of magnetic Weyl calculus implies that results depend on the properties of the magnetic field 
rather than those of the vector potential. 



1.2 Main results 

The main result of this work is Theorem II. II which gives an asymptotic two-parameter expansion of 
the product of two Hormander class symbols / e and g e S?™ 1 where each term involves the 
components of the magnetic field B, f and g and all its derivatives evaluated at the same point (x, £). 
Furthermore, two one-parameter expansions have been derived: for e <C 1, the expansion still has the 
same structure as the two-parameter expansion. In case e is not necessarily small, one can expand the 
product with respect to A <C 1 and the terms are given by oscillatory integrals. 

Theorem 1.1 (Asymptotic expansion of the magnetic Moyal product) Assume B is a magnetic field 
whose components are £% c £ 100 functions andf e 5" p mi as well as g e S^™ 2 . Then the magnetic Moyal product 
can be expanded asymptotically in e <C 1 and A <C 1: for every e<lwe can choose N = iV(e, e, A) e N 
such that 



nlxs=T,^ n ^Hff E , x g) M +h 



(1.8) 



n=0 k=0 



where the (n, k) term (/ttfxSOcn./fc) is in symbol class y™\ +m i O+*0p anc j we have explicit control over 
the remainder: R N as given by equation (I2.10D is numerically small and in the correct symbol class, 
5^ mi+m2 " (N+1)p , i. e. it is of order 0(e+) in the sense of Definitions^ The (n,k) term of the expansion, 



k \k 1 \---k n \' 

k o+2-, l= iJ k i =n 
ZU k i= k 



Y=X=Z 



is defined in terms of a family of differential operators S£j, j 
X {Y,Z) := \a(Y,Z) = \ (tj • z - y ■ Q 
1 



(1.9) 



X } {x,y,x). 



-ji S 3 - 

J m 1 ,...,mj_ 1 =l 



■ d x,„. B k iMy k z l -- 



c 



c = l 



s 

|a|+|/3|=j-l 



• ((1 - (-iy +1 )c - (1 - (-l) c )(; + 1)) y mi ■ ■■y mc _z I , 
C ha!p d x a dPB kl {x)y k z iy a zP. 



(1.10) 



6 



Here, we have glossed over the difficulty of agreeing up to which order we have to expand the product 
(e and A are independent), as we can no longer use well-known notation such as 0{e n ) or 0{X k ) without 
modification. We refer to Section |27l1 for details. 

Remark 1.2 With that in mind, we can give the first terms of the expansion concisely as 

d 

(f£,,8h, 0) = 4Z( v d ^ - v 

(=1 

d 

The second-order corrections contain at least two derivatives with respect to momentum; if we group by 
powers of e, then the decay properties are determined by (f j} B ,g)( 2 m> 

d 

U=i 

d 

If the magnetic field is constant, all terms containing derivatives of B vanish. 
For each order in e, the sum in A is finite and we immediately obtain the 

Corollary 1.3 (Expansion of |) B in e) If the assumptions of Theorem U.H are modified by taking X = \, 

then the e expansion of the product ft B eX g of two Hormander symbols can be recovered from the two- 
parameter expansion: the nth order term in e then reads 

k=0 

where the (J r (j^ %g\ n> k) ire taken from Theorem \l.l\ 

If e is not small, we can no longer expand the magnetic flux integral as a Taylor series in e, but have to 
accept it as-is. Only part of the twister can be expanded and the terms in the A expansion are obtained 
from replacing e lT '^ x,Y ' z ^ with 

e iT e ,,(.x,Y,Z) = e i§aCr.Z) e -iX^x,y,z) x e i|a(y,Z) ^ _ fryB^yrf + 0{X 2 )) 

in equation ( 11.61) . The integral formulas cannot be simplified any further unless the symbols which are 
multiplied with each other have a special structure (e. g. when they are polynomials in x or £). 

Theorem 1.4 (Expansion of |j B in A) Assume the magnetic field B has components of class 33 c £ co . Then 
for A <C 1 and e < 1, we can expand the A Weyl product of f e 5^ mi and g e 5" 1712 asymptotically in A such 
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that 



N 
Jc=0 

In particular, the zeroth-order term reduces to the non-magnetic Weyl product, (/ttf ja£)(o) = /tte£- We have 
explicit control over the remainder (equation (I2.12D ): if we expand the product up to Nth order in X, the 
remainder is of order 0{X N+l ) and in symbol class 5^i+™2-2(n+i)p 

The equivalence of the e — » X expansion to the X — > e expansion is obtained through explicit computation 
in Section 12.31 Agreeing on a remainder is somewhat tricky and necessitated the introduction of the 
concept of precision (Defintion |2.4D , because the numerical values of e and X vary independently. 

Theorem 1.5 Under the assumptions of Theorem U .1\ the magnetic Weyl product (P, of two symbols f e 
Sf™ 1 , g G 5^ m2 can be simultaneously expanded in e and X, i. e. the expansion is the same, regardless of 
whether we expand with respect to e first and then X or the other way around. 

1.3 Structure 

The derivation of our main results are found in Section [2} before we derive the main result, we need 
some prerequisites. First, the notions of two-parameter symbol classes and precision are introduced 
(Section [2TTT ). The properly adapted Wigner transform (Section [231 ) is necessary to show the equivalence 
of two versions of the product formula found in the literature (Section I2.3D . The one that is more 
amenable to an asymptotic expansion is used to derive the main result in Section \2A\ Lastly, we relate 
magnetic and non-magnetic quantization in Section [231 to be able to connect results derived via regular 
Weyl calculus to those where magnetic Weyl calculus has been used. 

As a simple, but non-trivial application, the semirelativistic Pauli equation is derived from the Dirac 
equation (Section [3]). It illustrates the versatility of the two-parameter expansion and gives insight into 
the origin of the corrections. We place emphasis on the mechanics of the computation. For the sake of 
brevity, the example is not presented in a mathematically rigorous manner, this is postponed to a future 
publication IIFL10II . 

In an attempt to clean up the presentation, we have moved some auxiliary technical lemmas and 
details of various straightforward, but tedious calculations to an Appendix. 

1.4 Acknowledgements 

This work was supported by a DAAD scholarship. I thank R. Littlejohn and R. Purice for their kind 
hospitality. I am very grateful for useful discussions, insights and references from M. Mantoiu, T. Miyao, 
G. Panati and H. Spohn. 

2 Asymptotic expansion in A and s 

This section will contain the proofs to my main results, namely the two-parameter expansion and some 
theorems which connect magnetic and non-magnetic Weyl calculus. Before we can attend to the asymp- 
totic expansion, we need some preliminaries: apart from assumptions on the magnetic field and some 
comments on the notation, we need to introduce the concept of precision as well as adapt the definition 
of the Wigner-Weyl transform. 

For simplicity, we will use Einstein's summation convention throughout this paper, i. e. repeated indices 
in a product are always summed over from 1 to d. We will always assume that the magnetic field 
satisfies the following assumptions unless explicitly stated otherwise. 
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Assumption 2.1 We assume that the components of the magnetic field B = dA and associated vector poten- 
tials A satisfy B kl e dS^ 00 ^ ,R) and A, e Sf p ~(R d ,R), respectively, for all l<k,l<d. 

Remark 2.2 If a magnetic field B satisfies the above assumption, it is always possible to choose a poly- 
nomially bounded vector potential, e. g. we may use the transversal gauge (equation (1B.3D ). It is also 
clear that if B and A satisfy this assumption, then so do the scaled field B e,A (x) := dA E ' A (x) = eAB(ex) 
and scaled potential A e - A (x) := AA(ex). 

In magnetic Weyl quantization, magnetic circulations and flux integrals play a very prominent role. We 
define the circulation of the one-form A along the line that connects x and y as 



r A ([x,y]) 



r 1 



A=(y-x). 

[*,y] 



dsA(x+s(y-x)). (2.1) 



The magnetic flux through the triangle with corners x, y and z (which we denote by (x,y,z)) is the 
(gauge-invariant) integral of the magnetic two-form, 



r B (x,y,z):= 



B. (2.2) 



Either we parametrize the triangle as in HIMP07II or we can choose a vector potential for B = dA and use 
the Stoke's Theorem to write r B (x,y,z) = r A ([x,y]) + r A ([y,z]) + r A ([z,x]). We will use the latter to 
derive the asymptotic expansion of the scaled flux integral 

T B (x,y,z) := T B ((x - |(y +z),x + f (y -z),x + §(y + *))) (2.3) 

in powers of e. 

2.1 Semiclassical symbols and precision 

The Hormander symbol classes are Frechet spaces whose topology can be defined by the usual family 
of seminorms 

l|/La : = SU P {K)- n+Hp \3^f{x,^\, fl ,ae< 

(x,E,)€T*R d 

One important notion is that of a semiclassical symbol llPST03bll . i. e. it is a symbol which admits an 
expansion in e and A which is in some sense uniform. 

Definition 2.3 (Semiclassical two-parameter symbol) A map f : [0, e ) x [0,A ) — > 5^ m , (e,A) — > 
f e,x is called semiclassical two-parameter symbol of order m with weight p e [0, 1], if there exists a 
sequence {/„, J n , fceNo , f n>k e ^ p m " (n+fc)p for all n, k e N , such that 

N 

(=0 rz+fe=f 

uniformly in the following sense: for each j e N there exists a constant C N m j > (independent of e and 
X) such that 



/ a -ty; £ n ^ k fn,k < C N>mJ max{e, A} 



(=0 n+Jc=! 

holds for all e e [0, e ) anc ^ ^ e [0, A ). 
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Since e and A vary independently, we also have to introduce a more sophisticated concept of precision. 
This is a technicality, but a definition is necessary to prove that expanding fif eX g first with respect to 
e and then A yields the same asymptotics as when the product is expanded with respect to A and then 
with respect to e (Theorem l2.13D . If there were only one small parameter, say e, then / — g = 0[e n ) for 
symbols /,ge implies two things: (i) the difference between / and g is 'numerically small' and (ii) 
we have associated a symbol class 5^ m ~ n P to the 'number' e". In case of two independent parameters, 
such a simple concept will not do and we have to introduce an association between a third number e <C 1 
and a certain symbol class. Although it seems artificial at first to introduce yet another small parameter, 
in physical applications, this is quite natural: say, we are interested in the dynamics generated by a two- 
parameter symbol H E ' A on times of order 0{ l /e), i. e. e~ l ' H '' . Then we need to include all terms in our 
expansion for which e n X k < e. Even if we choose e = e, for instance, we still cannot avoid this abstract 
definition as A is independent of e. 

Definition 2.4 (Precision ^(e+)) Let e < 1, A < 1. For e < I, we define n c ,k c ,N e N such that 

e n < +1 < e < e"s <e<X kc 



and N — N(e, A, e) := max{n c , k c }. We say that a finite resummation X„=o Sjc=o £n ^ k fn k °f a semiclassical 
symbol f E ' X e A5" p m is 0(e+~)-close, 

n=0 Jc=0 

ifff e ' X ~ Zto Zt ^ k fn,k e ^~^p and N E ,N X > N. 



2.2 Magnetic Wigner transform 

The Wigner transform plays a central role because it can be used to relate states (density operators) to 
pseudo-probability measures on phase space. We will need it to show the equivalence of two integral 
formulas for the magnetic Weyl product ^ x . 

Definition 2.5 (Magnetic Wigner transform) Let Lp,ip e 5^(R d ). The magnetic Wigner W A {f,ip) is 
defined as 

fr^ivMX) :=£ d {& a {y,w* x {-w)){-x). 

Lemma 2.6 The Wigner transform ifi"\((p,ip) with respect to g 5^(R d ) is given by 



2 - 

and maps 5"QR d ) x 5"QR d ) unitarily onto 5^(K 2d ). 

Proof Formally, the result follows from direct calculation. The second claim, Wf x (ip,ip) e 5^(M 2d ) 
follows from e -i^([*/s-y/2*/ e +y/2]) ^* (£ _ |) ^ (£ + 1) e x R d ) and the fact that the partial Fourier 
transformation is a unitary ony. □ 

Remark 2.7 The Wigner transform can be easily extended to a map from L 2 (R d x K d ) into L 2 (K 2d ) n 
< ig' 00 (R 2d ) where < g' 0O (R 2d ) is the space of continuous functions which decay at oo. For more details, see 
[|Fol89l Proposition 1.92], for example. 
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Let *i£^ lu (K ) be the space of smooth functions with uniform polynomial growth at infinity, i. e. for each 
/ e ^ lll (K 2d ) we can find mel,m>0, such that for all multiindices a,aeN^ there is a C aa > with 

\d2d*f(x,&\ <C aa {O m , V(x,£)eK 2d . 



Lemma 2.8 For y,^e 5"(R d ) and f 



3 polu 



(R 2d )C I ^ / (R 2d ) we have 



(27l) d 



Proof Since / e ^ p ~ lu (K 2d ) c ^"{R 2d ), it is in the magnetic Moyal algebra ^ B (K 2d ) defined in HMP04I 
Section VD.] and thus its quantization is a continuous operator 5^(IR d ) — > 5^(IR d ). Hence, the integral 
exists and we get the claim by direct computation. □ 

The Wigner transform also leads to a 'magnetic dequantization' - once we know the operator kernel, 
we can reconstruct the distribution. We do not strive for full generality here. In particular, unless the 
operator has special properties, we cannot conclude that / is in any Hormander class. More sophisticated 
techniques are needed, e. g. a Beals-type criterion IIIMP08II . 

Lemma 2.9 Assume B and A satisfy Assumption \2. 1 1 and Te3 (L 2 (IR d )) is a bounded linear operator 



whose operator kernel K T is in 



). Then the inverse magnetic quantization is given by 



dye 



2' 



+ !)• 



(2.4) 



This formula extends to operators with distributional kernels K T 
to continuous maps 5"(IR d ) — > 5^'(M d ). 



:5"( 



! ), i. e. the kernels associated 



Proof If T = Or/ A (/ r ) is the magnetic quantization of f T e 5^(M 2d ), then Wf x K T = f T e 5"(R 2d ) 
follows from direct calculation, using the explicit form of the Wigner transform, Lemma 12.61 Similarly, 
we confirm that T = 0^ eX (Wf x {K T )) holds and Wf x K T e ^(E 2d ) follows from K T e 5"(K d x K d ). 

If the kernel of T is a tempered distribution, then we can extend the formulas for Op^ x and 
to 5^'(M 2d ): Fourier transform, multiplication by a phase factor whose phase function is of tempered 
growth and a linear change of variables can all be extended to 5"'(IR 2d ) and thus it makes sense to write 
W^ X K T after a suitable reinterpretation. Then Wf x K T = f T e 5"'(R 2d ) is such that Op^fJV) = T : 
5^(M d ) — > y'(JSL d ). □ 



2.3 Equivalence of formulas for magnetic Weyl product 

It turns out that the integral formula for the product found in ||MP04l IIMP07II is not amenable to the 
derivation of an asymptotic expansion in e and A. Although an asymptotic expansion for e = 1 = A has 
been derived in HIMP07IL calculating each term has proven to be very tedious and it is not obvious how 
to collect terms of the same power in e and A. Thus, we will use an equivalent formula for the magnetic 
Weyl product. From this, we derive closed formulas for the {n, k) term by expanding the 'twister' of the 
convolution in the next section. 

Theorem 2.10 ( HMur99l[lMP07ID Assume the magnetic field B satisfies Assumption \2. 1 1 Then for two 
symbols f e S^™ 1 and g e 5" p m2 , the magnetic Weyl product f$ B eX g is in symbol class 5^ m i+ m 2 an d given by 
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the oscillatory integral 



(2ti) 



2d 



dY 



^a(* ~ fCK + + § (y - z),x + |(y + z)) (J^ 1 /)^) O^sK^) (2.5) 



(tte) 



2d 



dy 



dZ e-^ °<?-X'Z- x ) n B A (x - y + z, -x + y + z, x + y + z) f (?) g(Z) 



where Q B x (x,y,z) := e lXeT ^ x ^^ is the exponential of the magnetic flux through the triangle with corners 
x, y and z. 

Proof The Weyl product is defined implicitly by 

P ^(/)Op^(g)=:Op^(/tt^g) 

and its quantization maps 5^(R d ) to itself IIMP041 Proposition 21]. Combined with Theorem 12.91 this 
immediately yields 

where K op \ (/)o P A ( g ) is the kernel of Op EA (/)Op EA (g). Here, we have chosen a vector potential A which 
is associated to B that also satisfies Assumption 12.11 Although it is a priori not clear that there must 
exist a symbol f$ B x g, we will start with formal calculations and then use Corollary ID.3I to show that 
integral (I2.5D exists and is in the correct symbol class. 

Step 1: Rewrite in terms of Weyl system. Plugging in the definition of Op E A , we get 



(27l) 2d 

1 



(271) 



2d 



dy 



dy 



(2tt) 



2d 



• n B eX (Q e ,Q e + ey, Q s + ey + ez)W* x (Y + Z) 

dz[ [ dy C^; 1 / )m {K l s)(z ~ Y) e'^ [F - z) - 

■n B eX (Q e ,Q s + ey,Q e + ez) j W A A (Z). 



In order to find the kernel of this operator, we need to find the kernel for L E A (y,Z) := f2^ A (Q e ,Q E + 
ey,Q s + ez)W^ A (Z) which parametrically depends on y and Z = (z, Q. 

Step 2: Find the operator kernel for L eX (y, Z). Let <p e L 2 (R d ). Then we have 
(L E:A (y, Z)y>)(v) = n B x (ev, ev + ey, ev + ez) e-^+^H e-^Ctv.v+s]) ^ v + z ) 

due -ie(u-°/2)-n e -ar^[u-z,u]) n B x {eu - ez, eu + ey- ez, eu) 5 (u - (v + z)) y{u) 

duK LEX (y,Z;u,v)ip(u), 
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and we need to find Wf x K L ^ x {y,Z; •, -)(X), 



W? x K, s dy,Z;;-)(X). 



due-** e -'^«^-"/^"K ii£ , A (y,Z; f - \, *- + |) 
ia{X,Z) fiB ^ ( x _ | z , x _ £ Z + e y , x + | 2 ) =: La ( y) Z; x). 



Step 3: Magnetic composition law. Now we plug L eX {y,Z;X) back into the operator equation and 
obtain 



{2nf d 
1 



(27T) 



2<; 



dZ 



d7 



dZ e ^ + z) e ila{Y,z) ^ (x _ | (y + z)) x + | (y _ z)) x + s_ {y + z)) . 

•(^; 1 /)m(^; 1 g)(2). (2.6) 



This formula is the starting point for Miiller's and our derivation of the asymptotic expansion of the 
product. However, we can show the equivalence to the product formula obtained by two of the authors 
in HMP04II by writing out the symplectic Fourier transforms, 



RHS of 



(2ti) 



Ad 







c 


dY 


dY 


dZ 









d 2 e icr[X-Y,Y) £ ia(.X-Z,Z) & i e - a(J,Z) . 

n lx (* - §(y +*),* + fCy -*),* + §(y + z)) / (?) g(z). 



If one writes out the exponential prefactors explicitly, sorts all terms containing E, and 17 and then inte- 
grates over those variables, one obtains 



dF 



dZ e _i 7 x (y - g + X) y + z - X) -y + z + x) f (F) g(Z). 



Step 4: /jjg x g e t 5^ m i+ m 2_ ^he integral on the right-hand side of equation ( 12.61) satisfies the assump- 
tions of Lemma ID.3I with z = 1 = z' (keeping in mind that f2^ A satisfies the assumptions on G T < by 
Lemma [C2l) . Thus, the integral in equation d2.6D exists and is in symbol class 5* m i+ m 2. n 



2.4 Asymptotic expansion of the product 

To obtain an asymptotic expansion of the product, we adapt an idea by Folland to the present case 
[|Fol89l p 108 f.] : we expand the exponential of the twister 

00 00 

^ E E ^ E C nXa,a,b,pM y V*V 
n=0 k=0 

as a polynomial in y, 17, z and £ with coefficients C n j. a a b (3 e ^ < ^ 00 (K d ) that are bounded functions 
in x with bounded derivatives to all orders. Then we can rewrite equation ( 12.51) as a convolution of 
derivatives of / and g. Furthermore, we can show that there are always sufficiently many derivatives 
with respect to momenta so that each of the terms has the correct decay properties. 

The difficult part of the proof is to show the existence of certain oscillatory integrals. To clean up the 
presentation of the proof, we have moved these parts to Appendix [D] For simplicity, we also introduce 
the following nomenclature: 
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Definition 2.11 (Number of qs and ps) LetB e ^Sg" 00 ^, ¥™fiRtf x R^ d )) be afunction which can be 
decomposed into a finite sum of the form 

B(x,Y,Z) = b a abp(x,Y,Z)y a - n a z b ^ 

\a\+\b\=n 
\a\ + \P\=k 

where all b aab p smooth bounded functions that depend on the multiindices a,a,b,j3 e Njj. We then say 
that B has n qs (total number of factors in y and z) and k ps (total number of factors in 17 and Q. 

In the appendix we show how to convert qs into derivatives with respect to momentum and ps into 
derivatives with respect to position. Monomials of x and £ multiplied with the symplectic Fourier trans- 
form of a Schwarz function ip e Sf (R 2d ) can be written as the symplectic Fourier transform of derivatives 
of ip in £ and x: 

x a ^ a ip)(x) = ^((-id s r(id x ) a ip)(x) 

This manipulation can be made rigorous for symbols of Hormander class m with weight p . We see that 
derivatives with respect to momentum improve decay by p while those with respect to position do not 
alter the decay. In this sense, the decay properties of the integrals are determined by the number of qs 
and ps. 

Now we are in a position to prove the main result of this article, Theorem ll.il 

Proof (Theorem II. ID Let e <C 1. Then throughout the proof, we take JV — JV(e, e, A) e N to be as 
in the first part of Definition |2~4l i. e. e N+1 < e and A N+1 < e hold. We will show that ff eX g - 

Step 1: Formal expansion of the twister. We expand the exponential of the twister e I 2 t7(F ' z )e~ ,Ar * (xy,z) = 

e ir e A (x,y,z) up tQ t]le Ntil ternlj 

A/ 



e iT s ^x,m = £| L (r a(X) Y , z)) " + r n (x, y, z). 

n=0 

The remainder 



R N (x,Y,Z): 



1 rl 



dT(l-T) N 3 T N+1 e Tu| 



;N + 1 



r 1 



JV! 



(r EjA (x,7,z)) 



N+l 



dT(l-Tre lTJ ^ r ' ZJ (2.7) 



is treated in Step 3, right now we are only concerned with the first term. If we plug in the asymptotic 
expansion of the flux j B e derived in Lemma |RH up to JV'th order with JV' > JV, then we obtain 

(T EjA (x,F,Z))" = (fo-(F,Z) + AI:^ =le "> n ,(x,y, Z ) + AR iV ,[ r f](x,y, Z ))' , 

= 2(Y) { L 2 ^Y,Z) + XZ N n : =1 s n '^Ax,y,z)) n ~ l (XR N/ [r B E ]^,y,z)) 1 
1=0 ^ J 

=: (|a(y,Z) + A^' =ie "'^(x,y, Z )) n +R N , n [T EjA ](x,y,Z). (2.8) 
Again, we focus on the first term of the expansion and treat the remainder separately in Step 3: 
(|o-(y,Z)+A^' =ie "'^(x,y, Z )) n = 

= t £ e^xUa* - " ! , - (ja(r,# t n^(x, y ,,) 
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Now we define -Sf (Y, Z) := \<?{Y, Z) to clean up the presentation, include the sum over n again and sort 
by powers of e and A, 

N p 
n=0 

N i n n' n\ N ' 

*-An\ f- 1 fc !fc 1 !---fc Af /! \ f 

n=0 , u J. « , = i 

£j=o k i= n 



NN' n 



N' 



NN' n 



n=0 Jc=0 



n=0 ic=0 



Step 2: Existence of the (n, k) term. The properties of the (n, fc)th term of the product 



(.ff e ,,g\n,k)(X)- 



(271) 



2t! 



dY 



dZ e^ x - Y+ V p ( Xt Yy z) (J^gXZ) (2.9) 



can be deduced from the properties of ST n y. we proceed by showing that ST n k is a polynomial with x- 
dependent prefactors that contains n + k qs (powers of y and z) and at most n — kps (powers of 17 and 
0. .S? is the non-magnetic symplectic form and contains 1 q and 1 p. Hence, the fc th power of if 
contributes k Q qs and an equal amount of ps. By Lemma IBTTl the magnetic terms if,-, j > 1, contribute 
j + 1 qs and no ps. In this sense, magnetic terms improve decay. By conditions imposed on the indices 
appearing in the definition of SF n k , we deduce there are 

N' N' N' 

fc + Yfi + 1)fc j = fc ° + Ti^j +Y i k J =n+k 

j=l ;=1 j=l 

qs and fc ps. As < fc < n — k, Lemma [LT2l implies the existence of integral J2.9I ) and that it belongs to 
the correct symbol class, namely (/j)^ x g\ n ,k) G ^i+^-^+^p. 

Step 3: Existence of remainders. There are two remainders we need to control, equations ( 12.7D and 
02.8D : the first one stems from the Taylor expansion of the exponential, the second one has its origins in 
the expansion of the magnetic flux, 



R= (x, Y, Z) := R N (x, Y,Z) + ^ -R H , n [T s>x ] (x, Y, Z). 



The remainder of the product is obtained after integration, 



R N (X) :-- 



(271) 



2d 



dy 



dZe i^,Y+z) R ^ Y>z) (Jf; 1 g)(Z). (2.10) 



We have to show that (i) the integral exists, (ii) it is in the correct symbol class and (iii) it is of the right 
order in e and A. Points (i) and (ii) are the content of Lemma ID . 3 1 and we have to show that each of the 
two contributions to the remainder satisfies the assumptions. 

The first contribution to R N stems from the Taylor expansion of the exponential, 



(2*0 



2d 



dY 



r 1 



dZe 



ia(X,Y+Z) 



dT(1 _ T) N a N + l e ,u| (Jr - 1/)(y)(J ,-l gXz): 



r 1 



(2ti) 



2d 



dT(l-T) A 



dY 



,x(x,Y,Z) ' 
,-N+l 



dZe^ x - Y+z ^—{T e ^x,Y,Z)) N+1 e- l ^r>,y^. 



,iTf<r(y,z) 



c^/xnc^gxz). 
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The first factor, (T E Y,Z)) N+1 , can be expanded in powers of a(Y,Z) and j B (x,y,z): 

N+l , s 

(r E , A U,y,z)) N+1 = y ~i j a ' (^(^r^'Mx^y 

As e N+1 < e holds by definition of JV, the first term of the remainder is of the correct order. The decay 
properties are dominated by (cf(Y,Z)) n+1 with JV + 1 ps and N + 1 qs. All other terms contribute less 
than JV + 1 ps and more than JV + 1 qs since by Lemma iBTTl y E is of order e and contributes 2 qs and no 
ps. Furthermore, Lemma ICTl gives polynomial bounds of derivatives of y E : 

\d x a r B e ix,y,z)\<C a ((y) + (z)) lal 
A similar bound holds for the exponential of the flux ("Corollary IC.2D : 

\d«e- lX ^ x ^Xx,y,z)\ < C a (y) |a| (z) |Q| VaeNj 

Altogether, (T e ^x,Y,Z)) N+1 e -^^r B A^,y,^ satisfies the conditions on G T , in Lemma ILT31 (with t = t') 
which implies 



(271) 



2d 



r 1 r r ,-jv+i 

dZ e ia *- T+z > — (T E _ A (x, Y, z)) N+1 e- M ri(*,y*\ 



dT(l-T) 





N 



dY 



■ e 2 



exists as an oscillatory integral and belongs to symbol class 5^ m i +m 2 (w+i)p_ 
The second contribution which involves 

R N 'n [T e> J (x, Y, Z) = J] " (f a(Y, Z) + A£*' = i e n ' ££ n , (x, y, zj) ( XR N , [ Y E ] {x, y, z)) 1 
i=i ^ ' 

can be estimated analogously: by Lemma [B7TI R N >[y B ] is of order 0{e n +1 ) (the largest prefactor is 
e N +1 A < e) and contains JV' + 2 qs. So the terms in the above sum contain at least JV' + 1 > JV + 1 more 
qs than ps and another application of Lemma ID.3I (with t = 0) implies that the second contribution to 



R N exists as an oscillator integral and is of symbol class c 5^" 1 +m 2 -(jv'+i) p c 5 ^i+m 2 -(w+i) P 
Altogether, we conclude that R N exists pointwise, is of symbol class 5^ m i+ m 2-(w+i)p as i on g as JV' > JV 
and hence ff eX g - X^ =0 £Lo e n X k {ft e>x S\n,k) = &{e+). This concludes the proof. n 

If we do not have a separation of spatial scales, i. e. s = 1, but weak coupling to the magnetic field, 
we can still expand the product fl^ A as a power series in A. This is also the starting point of the A-e 
expansion which coincides with the e-A expansion. 



Theorem 2.12 Assume the magnetic field B satisfies Assumption \2. 1 \ then for A <C 1 and s < 1, we can 
expand the A Weyl product of f e 5f™ 1 and g e 5^ m2 asymptotically in A such that 

nix* - |>*cf e 

fc=0 

where {f^ e}i s\k) e 5^jn 1 +m 2 -2kp ^ gj ven fry equation ( 12. 1 ID . In particular, the zeroth-order term reduces 
to the non-magnetic WeyZ product, (/jF jg)(o) = /tleg- We ^ ave explicit control over the remainder (equa- 
tion ( 12.12D ): i/ we expand the product up to Nth order in A, the remainder is of order 0{X N+l ) and in 
symbol class i5 ^ 1 + m2 -2(N+i) P ^ 
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Proof Assume we want to expand up to Nth order in A. We will show f$ B , g — 2k=o ^ k (f x&\k) = 
6{X N+l ) and that the difference is in 5 ^ 1 +m 2 -2(N+i) P _ 

Step 1: Expansion of exponential flux. If e is not necessarily small, we cannot expand the magnetic 
flux integral j B in powers of e anymore. However, we will keep e as a bookkeeping device. Expanding 
the exponential of the magnetic flux, we get 

e iT t> Xx,Y,Z) _ e qa(Y,Z) e -iXf e {x,y,z) 



i a(Y ' Z) (iLo? if. ix,y,z)) k +R N {x,y,z)). 



The remainder is of order A + and has 2(JV + 1) qs, 



R N [x,y,z)= —(-iXy B {x,y,z)) 



N + l 



This can be seen more readily once we define — eB^.[x,y,z)y l Zj := j B {x,y,z) to emphasize that y B 
contains e as a prefactor and 2 qs. Using the antisymmetry of B (j , there is a simple explicit expression 
for B*. (see proof of Lemma [BTTb : 



/•+V2 



%(x,y,z) = - 



r 1 



dt 



1/2 



dss [B l} (x + es(ty- */2)) +B lj (x + es(y/2 + tz))] = 0{l) 



Step 2: Existence of fcth-order term. Then the expansion can be rewritten so that we can separate off 
factors of y , z and e. The fcth order term contains 2k qs and no ps, 



(-0* 

fc! 



;k k 



f s ix,y,z)y 



fc! 



l^ iJm (x,y,z)y, m z jin 



By Lemma IL\3l (with z = 1 = z') the fcth order term of the product 

J<- r r r ;k k 



(2*0 



2d 



dY 



•G^Vxno^gXz) 



(2?r) 



2d 



dy 



dZ e to »' y+ ^ e 1 ^ - Y\Bl jm (x,y,z) ■ 

^ " m=l ' 



• •■■\/)Kn(^; 1 (%. -9 t g)){z) (2.ii) 



exists and is of symbol class ^^i+ m 2 2k P. 

Step 3: Existence of remainder. The remainder is of order A N+1 and has 2(JV + 1) qs. It contains e N+1 
as a prefactor as well which will be of importance in the proof of the next theorem. By Lemma [CTl and 
Corollarv lC.21 the integral in R N over the exponential of the magnetic flux is bounded and its derivatives 
can be bounded polynomially in y and z, 



„N+1 



r 1 



N+l 



■{B e l Xx,y,z)y l z ] ) 



dz'{l-z'f e-'' AT '^(*,y,z). 



This means R N satisfies the conditions on G z > in Lemma HOI (with z = 1) and we conclude that 
1 



R N (X) r- 



dY 



dZe 



ia(X,Y+Z) J^cr{Y,Z) 



e [ 5 CT(y ' z) R N (x,y,z)(^ ff - 1 /)(7)(^ (T - 1 gXZ) (2.12) 



exists and is in symbol class 5^" 1 +m 2 -2(iv+i)p_ 
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The next statement is central to this paper, because it tells us we can speak of the two-parameter expan- 
sion of the product. 

Theorem 2.13 Assume that the magnetic field B satisfies Assumption ^. l\ and e <C 1 in addition to X <C 1. 
Then we can expand each term of the X expansion o//(j^ A g in e, f e S^ mi , g e 5*"" 2 , and obtain the same 



as in Theorem \l.l\ Hence we can speak of the two-parameter expansion of the product 



Proof Step 1: Precision of expansion. Assume we have expanded the magnetic product fl^ up to 
JV th power in X. Then for the remainder of the proof, we fix N = N(e, s, X) 6 N as in Definition [Z4] for 

e = X N ° . 

Step 2: Equality of (n,fc) terms of expansion. Now to the expansion itself. The two terms we need 
to expand are the non-magnetic twister e l z a ^ Y ' z ^ and the fcth power of the magnetic flux integral j B e in 
£<1: we choose N',N" > N and write the fcth order of the X expansion as 



(2*0 



2d 



dY 



(2k) 



2d 



dY 
(-if 



dZ e-C^+iO ^ QS n L(l a( Y,Z)) n +R N ,l<T](Y,zj)- 



The remainders are given explicitly in Step 3, equations A2.13D and H2.14D . The (n, fc) terms of the 
expansion originate from the first of these terms, i. e. we need to look at 



N' 



nlfcj ! • • • fc N "! 



to obtain the (n, fc) term of this expansion. The remaining three terms define the remainder which will 
be treated in the last step. We define if (7,Z) := ^<j(Y,Z), fc := n and recognize the result from 
Theoremll.il the terms match: 



N'N 

I 

n=k 



N" 



fc !fc 1 !---fc N » 



;=i 



Obviously, the arguments made in the proof of Theorem fTTTI can be applied here as well, and we conclude 
that the (n, fc) term exists and is in the correct symbol class, 5^ m i+ m 2-(n+'0p_ 

Step 3: Existence of remainders. The remainders of the expansions of e 1 2 <7(y,z:) and (y B (x,y,z)) k , 



■N' + l 

iV[a](Y,Z) = ^' +1 — (±a(7,Z)) N ' +1 



r 1 



dT(l-T) W V'i T<7(r ' Z) 



(2.13) 



and 



R N » k [XR](x,y,z) = Yl(j) (Sji^Cx.y,*))*" 1 (^[ r f](x,y, Z ))' 



(2.14) 
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with Rjv"[yf](x, y,z) as in Lemma [570 lead to three terms in the total remainder: 

Rl N , N „ k (x, Y,Z)=R N ,[cr](Y,Z)( (x, y, z)) k + R N „ k [S£R] (x, y, z)) + 

+ (S!!' : loe n £(^a2))")R Nn [^R](x,y, Z ) 

Going through the motions of the proof to Theorem ll.il we count ps and qs, and then apply Lemma IL\3l 
The first remainder, R N >[a][Y,Z), is of order e N+1 < e in e and contributes N' + 1 qs and ps. By 
Lemma [BTT1 R N » [yf ] contributes at least N" + 2 qs and all prefactors are less than or equal to e N +1 < e. 
Thus the terms in R N « k [J£R] contain at least N" + 2 qs (for all k < N) and prefactors that are at most 
e N +1 < e. Hence, the total remainder exists as an oscillatory integral, is 0[e+) small and in symbol 

Class 5"ni+m 2 -(N+l)p_ q 

Remark 2.14 The asymptotic expansion of ^ x can be immediately extended to an expansion of products 
of semiclassical two-parameter symbols (see Definition l2.3t . 



2.5 Relation between magnetic and ordinary Weyl calculus 

In a previous work [IMP07], Iftimie et al have investigated the relation between magnetic Weyl quanti- 
zation and regular Weyl quantization combined with minimal substitution, the 'usual' recipe to couple a 
quantum system to a magnetic field. However, since there were no small parameters e and A, we have 
to revisit their statements and adapt them to the present case. 

Let us define itf (X) := £; — AA(x) as coordinate transformation which relates momentum and kinetic 
momentum. With a little abuse of notation, we will also use / o := /(x,#^(X)) to transform 

functions. In general, Op^ ,(/) ^ Op e (/ o since the latter is not manifestly covariant. However, we 
would like to be able to compare results obtained with magnetic Weyl calculus to those obtained with 
usual Weyl calculus and minimal substitution. To show how the two calculi are connected, we need to 
make slightly stronger assumptions on the magnetic vector potential. This may appear contrary to the 
spirit of the rest of the paper where it has been emphasized that restrictions should be placed on the 
magnetic fie Id. The necessity arises, because usual, non-magnetic Weyl calculus is used in this section. 

Assumption 2.15 We assume that the magnetic field B is such that we can find a vector potential A whose 
components satisfy 

\d x a Mx)\ <C a , VI < l < d, \a\ >l,o6 Nq. 



In particular, this implies that the magnetic field B = dA satisfies Assumption 12.11 i. e. its components 
are gS'tf 00 functions. It is conceptually useful to introduce the line integral 



r 1 

rAr 



dsA(x+s(y -x)) (2.15) 

o 

which is related to the circulation r A ([x,y]) = (y — x) ■ T A (x,y); similarly, T A ([x,y]) =: (y — x) • 
T A (x,y) defines the scaled line integral. This allows us to rewrite the integral kernel of a magnetic 

/ A (/)for/ e5 ; m 



T A (x,y) := 



pseudodifferential operator Op A ,(/) for / e5"" as 



Kf,e,*.(x,y) 



dr? e-^f (f (x + y), r, - AT A (x, y)) . (2.16) 



If we had used minimal substitution instead, then we would have to replace the line integral r A (x,y) by 
its mid-point value A(|(x + y)). 
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Theorem 2.16 f HIMP07I0 Assume the magnetic field satisfies Assumption \2.15\ Then for any f <E 5" p m 
there exists a unique g e 5^ m such that Opy^f/) = Op E (g °'0 A )- g can be expressed as an asymptotic series 
g x IZ o XL: Sn,k, ^ere g nM e 



/or aZZ n > 1, and 



2 A fc g n , fc (x, O = E -"^i (id y ) a [pif (x, I - XY\x + ly, X -Zy) + AA(x))) 

fc=l a|=n a " 



y=0 



(2.17) 



OnZy terms with even powers of e contribute, i. e. g n k = for all n e 2N + 1, 1 < k < n. In particular we 
have g ao = /, g 1>0 = 0, gi A = and f - g e 5*^-^. 

Remark 2.17 The reason that only even powers of e contribute can be traced back to the symmetry of 
r^(x, y) = +r^(j,x). Note that this is consistent with what was said in the introduction, r A ([x,y]) = 
(j — x) ■ r A ([x,y]) is indeed odd. 

Proof The proof is virtually identical to the proof of Proposition 6.7 in HIMP07II ; we will only specialize 
the formal part to the present case, the rigorous justification can be found in the reference. 

For a symbol / e 5^ m , the integral kernel of its magnetic quantization is given by equation ( I2.16D . 
On the other hand, it is clear how to invert Or/ A for A = 0, A = 0: we apply the non-magnetic Wigner 
transform W e := ^ff£ to the magnetic integral kernel of Op^ A (/): 



dye-^ Wf + f, f -f) 



dr, e^ f (x, r, + £ - Ar A (x + § y, x - f y)) 



Since we have a separation of scales, we can expand r A (x + |y,x — |y) in powers of e up to some 
even N. We will find that only even powers of e survive - which immediately explains the absence of the 
first-order correction, 



r A (x + , x , ) 



ds 



2> v 

n=0 al=n 



n=U v ^ a =2n 



3 x °A(x)y a +R N (s,x,y) 



dsR N (s,x,y). 



■1/2 



The remainder is bounded since it is the integral of a function over the compact set [— 1/2, + 1 / 2 ] x 
[0, 1] . In any event, The exact value will not matter if we choose JV large enough as we set y = in the 
end. 

A Taylor expansion of / (x, 17 + E, — Ar A (x + ~y,x — §y)) around 17 — Ar A and some elementary 
integral manipulations formally yield for the nth term of the expansion 



1 y] A fc g njk (x, ? - AAfx)) = ~, (idyY (x, £ - Ar A (x + § y,x - § y))) 



/t=i 



|a|=n 



y=0 



(2.18) 



where we substitute the expansion above for T A . Each derivative in y will give one factor of s, i. e. we 
will have n altogether. On the other hand, we have at least 1 and at most n factors of A. Only even 
powers in e contribute, because the expansion of r A (x + |y,x — |y) contains only even powers of y. 
Furthermore, all terms in this sum are bounded functions in x, because derivatives of A are bounded by 
assumption. 

To show that g n k is in symbol class 5^ m -( n+i 0p j we need to have a closer look at equation (12.18D : the 
only possibility to get k factors of A is to derive dff (x, £ — Af A (x + | y, x — ~ y)) k times with respect 



20 



to y. Each of these y derivatives becomes an additional derivative of d?f with respect to momentum. 
Hence, there is a total of \a\ + k = n + k derivatives with respect to £. 

The rigorous justification that these integrals exist can be found in HIMP071 Proposition 6.7]. n 
Remark 2.18 If we are interested in a one-parameter expansion in e only, then 

g n {X) := e- n ^ W (^7 (*, ? - Ar A (x + | y,x - |y) + AA(x))) | _ q 

|a|=n 

gives the nth order correction in e. 

Proposition 2.19 ([IMP07]) The converse statement also holds: if the magnetic field satisfies Assump- 
tion \2.15\ then for each g e 5" p m there exists a unique f e such that Op E (g o = Op^ A (/), 

/ x 2^=o 2fc=i en ^ k fn,k> fn,k e t 9^ m -("+ /c )P ; can foe expressed as a formal power series in e where the 
nth term is given by 

J X k f nM {x, £) = e-"^ (i3 y ) a (S«f) (x, £ + XT%x - y/2, x + r/2) - AA(x)) | y=0 (2.19) 

fc=l |a|=n a ' 

In particular we have / = g, / 10 = 0, = and g — / <E ^™- 3 P m 

Proof This proof works along the same lines: one magnetically Wigner- transforms the kernel of the 
operator Op E (/ o #^), we refer to [TMPOZ, Proposition 6.9] for details. □ 



3 Application to the Dirac equation 

To demonstrate the advantages of magnetic Weyl calculus, we will apply it to a simple, yet interesting 
problem: the semirelativistic limit of the Dirac equation. This is a well-studied problem IIFW50l|Tha921 
|Ynd961 ICor831 [Cor04ll , but we believe our derivation sheds a new light on origin of corrections. To keep 
this section readable and put emphasis on the computational aspects, we will dispense with mathemat- 
ical rigor. Making these statements exact and putting them into context with previous works will be 
postponed to a future publication IIFL10II . 

The dynamics of a relativistic spin-!/2 particle with mass m subjected to an electromagnetic field is 
described by the Dirac hamiltonian, 

id t V = (c 2 mf3+c (-iVj • a - eA(x) • a + eV(x) j *, * e L 2 (M 3 , C 4 ). 

The hamiltonian consists of operator-valued matrices: a ; , j = 1,2,3, has the ;th Pauli matrix as entries 
in the offdiagonal, [3 is the diagonal matrix with entries 1, 1, —1 and —1, namely, 

P { -id c2 J • 

As is customary, we have used shorthand notation for £, ■ a := Xj=i a j- If we assume that the compo- 
nents of the magnetic field B = dA and an associated vector potential A satisfy Assumption 12 . 1 1 and that 
V e 33 S^QR 3 ), then Proposition 1.1 from HHNW89II as well as the fact that multiplication by V defines 
a bounded operator on L 2 (R 3 ,C 4 ), the Dirac hamiltonian defines an essentially selfajoint operator on 

^ C °°(1R 3 ,C 4 ). 

If we rescale the energy by 1 /c 2 for convenience and absorb the charge into the definition of the 
potentials, we see that there are two natural ways to write the Dirac hamiltonian, namely 

H D = m/3 + i (-iV x - ±A(Q)) • a + ^V(Q) 
= m/3 + (-i±V x - ^A(Q)) • a + ± V(Q) 
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where Q := x is the position operator. The first way of writing suggests to use 

p? := -iV x - ~A[Q) 

as kinetic momentum operator, the second definition, 

P A c := -flv, - ^A(Q) = Ip* (3.1) 

absorbs an additional factor of l jc. This seems nothing more than an algebraic trivially, but is is this 
choice of building block operators which contains the physics. The first corresponds to the non-relativistic 
scaling where momenta are very small and the 1 /c — » limit leads to the non-relativistic limit. In IIFL10II 
we derive the Pauli equation including fourth-order corrections with P^ 2 as kinetic energy operator 
in this scaling. This case is computationally more involved, than the second, the semirelativistic scaling, 
which we will discuss now. The limit l /c — » will lead to the semirelativistic quantization where the 
kinetic energy operator is the magnetic quantization of \/m 2 + E 2 . The rescaled Dirac hamiltonian can 
be written as 

Hd=H (P?)+^H 2 (Q) (3.2) 

where 

H 2 {x) := V(x). 

Or to put another way, we can write H D = Op^(H D ) as the magnetic quantization of the symbol H D := 
H Q + -^H 2 with respect to the pair of observables (Q, P^). The attentive reader will notice that we have 
defined the magnetic quantization of matrix-valued symbols; we account for this by tensoring the Weyl 
system with the identity matrix id C 4, 

dXJf- 1 (H D )e iCT ^ (Q ^®id C 4. 

To put things into perspective, let us mention some works pertaining to the current approach: most of 
the previous treatments (e. g. ||Hun75| [GNP89I |Tha92l |Ynd96|| ) are dealing with the derivation of the 
non-relativistic limit where in our language the usual Pauli hamiltonian is the effective hamiltonian. Note 
that some of the authors use a different scaling to derive the non-relativistic limit: the prefactor of A is 
taken to be c and not 1 in the original Dirac equation. In that case, the magnetic and electric field scale 
differently. 

The first to apply pseudodifferential techniques to the Dirac equation in order to obtain the semirela- 
tivistic limit was Cordes IICor83ll who ordered the terms of the diagonalized hamiltonian by decay and 
not by powers of a small parameter. Physically, this is not satisfactory, because the prefactor decides 
which effects are and which are not measurable. Furthermore, it is not clear how to extend his ideas to 
allow for a non-relativistic limit. Brummelhuis and Nourrigat HBN99II also consider the problem with the 
help of pseudodifferential theory they derive a power expansion in the semiclassical parameter ft for ap- 
proximate projections onto electronic and positronic state as well as the effective hamiltonian. They give 
the first-order correction to the effective hamiltonian explicitly and then continue with a semiclassical 
limit via an Egorov theorem. 

We adapt a technique developed by Panati, Spohn and Teufel, space-adiabatic perturbation theory 
IIPST03bll . invented to cope with more general problems and due to its simplicity, the Dirac equation is 
one of the first systems this technique has been applied to IITeu03l Section 4.1]. Their derivation rests on 
a very different physical mechanism, though: in their approach, the adiabatic decoupling is due to a slow 
variation of the electromagnetic potentials and in the limit e — » 0, the particle needs to travel farther 
and farther to see appreciable differences in the potentials. This limit is in fact equivalent to the case 



Op?(H D ) := 



(2n) 3 
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considered by Brummelhuis and Nourrigat HBN99II (see Appendix©. We, on the other hand, assume 
that the relativistic energy is small compared to the rest energy mc 2 of the particle and consequently, 
the particle's velocity is small, v o/ c < 1- While this sounds very similar to the point of view of Panati, 
Spohn and Teufel, the fields scale differently: the electromagnetic fields associated to the slowly-varying 
potentials A(ex) and V(ex) are of the order 0{e) while the electromagnetic fields in our case are of 
order as we will see. Furthermore, we do not need to introduce another scale (the slow variation 

of the potentials) in addition to the energy scale given by l jc. Our result holds almost in the entire range 
of validity of the Dirac equation: if the total energy of the particle approaches the pair creation threshold 
2mc 2 , we no longer expect the Dirac equation to give an accurate description of the physics anyway. 

Before we proceed, we give the first few terms of the asymptotic expansion in 1 /c of the magnetic Weyl 
product $ . 

3.1 Asymptotic expansion of ^ 

If we compare equation HI. ID with the definition of P^, equation J3.1D . we see that e = l /c and A = 1 /c 2 . 
According to a simple modification of Theorem IA.3I we can write the expansion of in terms of the 
two-parameter expansion of (J^ , . The first few terms of / j)^ g (with / and g being suitable matrix-valued 
functions, e. g. matrix-valued Hormander class symbols) are 

(ff e g\o)=fg, (3.3) 
iff c g\i) = - L 2 {f,g}, 

(ff c gh) = - l MV Y ,V z )ff{Y)g{Z)\ Y=x=z , 

(ff c gh) = L M^y, Vz)Yf(Y)g(Z% =x=z + |B iy (x)%/(X)3^fX). 

While this seems very complicated, we will often need the product of two symbols which are functions 
of momentum only, / = /(£), g = hi that case, only purely magnetic terms (i. e. fc = in 

Theorem ll.lD contribute, 

ff c g=fg + ^ L 2 *iAj %g + wo =fg+ wo. 0.4) 

3.2 Semirelativistic limit as adiabatic limit 

The technique of choice, a modified version of space-adiabatic perturbation theory [|PST03bl IPST03al 
ITeu03ll that uses magnetic Weyl calculus, rests on the interpretation of the semirelativistic limit l /c —* 
as an adiabatic limit. This means, the Dirac hamiltonian has three characteristic features all adiabatic 
systems share, the so-called adiabatic trinity: 

(i) A distinction between slow and fast degrees of freedom, i. e. a decomposition of the original Hilbert 
space the hamiltonian acts on into #6 = ^iow®-^fast- Here, the fast Hilbert space is spanned by the 
electronic and the positronic state, Jtf ast = C 2 . The slow Hilbert space is that of a non-relativistic 
spin- 1/2 particle, J«f slow = L 2 (K 3 ,C 2 ). 

(ii) A small, dimensionless parameter that quantifies the separation of scales. If v is a typical velocity 
of the particle, we expect that no electron-positron pairs are created as long as v o/ c ^ 1- However, 
for notational simplicity, we use 1 /c as small parameter. 

(iii) A relevant part of the spectrum of the unperturbed operator, separated by a gap from the remainder. If 
we consider the field-free case, then H (— l /c V x ) fibers via the Fourier transform and the spectrum 

of each fiber hamiltonian is given by spec(H (£)) = {±^/m 2 + ^ 2 }. We are interested in the 

electronic subspace - which is separated by a gap (of size 2ym 2 + ^ 2 > 2m) from the positronic 
subspace. This ensures that even in the perturbed case, transitions from one band to the other are 
exponentially suppressed. 
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In a commutative diagram, the unperturbed situation looks as follows: 



-itffoC-tycV*) 

n 



itM(.-'/c V x )fl 

n 



L 2 (K 3 , C 4 ) "° ( i/cV,) : L 2 (K 3 , C 2 ) ® C 2 



(3.5) 



rcoHAV x )(L 2 (R 3 ,C 4 )) 



^L 2 (K 3 ,C 2 ) 

u 



With a little abuse of notation, we will interpret all of these spaces as (subspaces of) L 2 (R 3 ,C 4 ) = 
L 2 (R 3 ) <g> C 4 when convenient; operators acting on this space can be thought of as (operator- valued) 
4x4 matrices or, if we are on the right-hand side of the diagram, as 2 x 2 matrices whose entries are 
itself (operator-valued) 2x2 matrices. The former identification is used during calculations, but the 
latter is conceptually useful. 

The objects in this diagram can be found in every text book on relativistic quantum mechanics (e. g. ||Tha92[ 
|Ynd96|| ): tz is the projection onto the electronic subspace, 



7 r (?) = i('id C 4 + i ^H (?) 



u is the matrix- valued function that diagonalizes H , 



h := u H u * = ^m 2 + ^l5 =: E/3, 
and 'intertwines' n with the reference projection, 



i ;l u 



L ref " 



id C 2 




where 



y/2E(E + m) 



((E + m)idc4-(?-a)j8). 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



The quantization of n re{ is n ref = id L 2 (K 3 C 2) ® 7i ref projects out the positronic degrees of freedom in 
diagram J3.5D . If we are interested in the electron's dynamics only, we can describe it by an effective 
hamiltonian, the quantization of 



(3.10) 



KffO := ^ref h 7I ref = E id C 2 = \j W 2 + E, 2 id C 2, 

in the following sense: 

(e- ltH » H/<7 - ) - u *HA Vje-' 1J (-^»'u (-./c Vj) *„(-</« VJ = 

Hence, we are able to relate the dynamics in the upper-left corner of diagram 03.5D with the reduced, 
effective dynamics in the lower-right corner. This reduction is possible as H (— '/c V x ) and 7i (— ; /c V x ) 
commute, [H (— i/c V x ), n (— i/cV,)] = 0. Hence, the electronic subspace is invariant under the unper- 
turbed dynamics. 
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If we switch on the electromagnetic perturbation, this is no longer true, the commutator of H D and 
Op^(7i ) = 7i (P^) is °f order The immediate question is whether we can generalize diagram 

J3.5D through some generalized projection IT and generalized unitary U c such that 



D 

L 2 (K 3 ,C 4 ) 



n c (L 2 (K 3 ,C 4 )) 



C) 

L 2 (R 3 , C 2 ) ® C 2 



^L 2 (K 3 ,C 2 ) 



(3.11) 



holds. If these objects exist, we require them to be an orthogonal projection and a unitary which commute 
with the full, perturbed Hamiltonian H D and block-diagonalize it, i. e. 



n" = IT, n 

7 C* TjC 



rT 



u c *u c = id L 2 (M , u c u c 



id. 



[o P ^(H D ),n c ] =o 
u c rr u c * = n rPf = id r2r „ 3) ® 7i rpf . 



Because of the last property U c , is called intertwiner. For suitable potentials V, we can translate these 
equations (up to A5^~°°) into equations of semiclassical symbols. So if there exist 7r c e A5^° and u c e A5^° 
such that IT = Op^(n; c ) + ^oCV^fl and U c = Op^(u c ) + ^bOA "), then the corresponding symbols must 
satisfy 



[H D ,7i c ] 8B = ^(Vc») 



id C 4 + <?(i/c»), iz c tty * = id C 4 + <7(i/c-) 



u nn nil 



(3.12) 
(3.13) 



where the Moyal commutator is defined by \_H D , tt c ] u B := H D ^n c — tz c ^H d . 

If we incorporate magnetic Weyl calculus into space-adiabatic perturbation theory IIPST03bllTeu03ll . 
we obtain an explicit resummation for these symbols as well as formulas to correct Op^(n: c ) and Op^(u c ) 
on the order ^(Vc") t0 get a tme projection and a true unitary in the operator sense. Then it is natural 
to assume that the principal symbols (the zeroth-order term) of the expansion of if and u c have to be n Q 
and u - the symbols of unitary and projection associated to the unperturbed hamiltonian. Starting from 
the unperturbed objects, Panati, Spohn and Teufel have found recursion relations which give corrections 
to u and n order-by-order in 1 /c which turn out to be independent of the specific Weyl calculus used. 

The generator of the dynamics in the lower-right corner of diagram (I3.11D is the upper-left 2x2 
submatrix of the diagonalized hamiltonian 



i. e. the effective hamiltonian 



Kff ■= n rel hTi re( = 7i ref u c f_ H D fu c * 7I ref . 



(3.14) 



(3.15) 



There are technical and conceptual reasons for this specific choice that go beyond the scope of this text, 
we refer the interested reader to IITeu031 Section 3.3] for details. The magnetic quantization of h efS 
generates effective dynamics which approximate the full dynamics for electronic states, 

r e -itotf{H») _ op^( u c )* e - lt0 ^ h ^ Op^(u c )) Opf(7i c ) = C {\t\ i/c»). 



: We say that two c-dependent bounded operators A and B on a Hilbert space J^f satisfy A = B + ffo( 1 / c °°) if f° r eacn n s N there 
exists a constant C n such that HA-BH^r^) < C n \ 
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3.3 Effective hamiltonian 

In the present case, equation ( 13.41 ) implies that the first correction to n and u Q is of third order in l /c\ 
TTottf ^0 - n = 0(VcO [H D , n ] {f = ( Ye 3 ) 

U 1> = id C 4 + 0^), U f c U * = id C 4 + 0(l/c 3 ) U j)> t)f Uq* = 7I ref + ^(Vc 3 ) 

From these equations, we conclude that n Q and u are an approximate Moyal projection and Moyal 
unitary, respectively, i. e. n c = n Q + and u c =u + ^(Vc 3 ). 

To compute the terms in the expansion of h eif up to third order in 1 /c, we need to obtain the diagonal- 
ized hamiltonian symbol h := u c ^H D ^u c * up to second order first. As expected, the leading-order term 
is the relativistic kinetic energy, 

h = (u f c H f c u *) m =u H u * = Ep. 

If ft ~ Xi^Lo jzhn i s the asymptotic expansion of the diagonalized hamiltonian, then we can determine h n 
recursively from h^u c = u c ^H D + ff^/c 00 ): 

7^jf" c ) ( o) = ^„"o + ^( 1 A" +1 ) 

= " c ttfH D - (JZl 7 h)fy + 

This simplifies calculations considerably. Starting from this equation, we arrive at the following formulas 
for h 1 and h 2 : 

hi = (u H x +u 1 H - h Ul + Cuottf H \i) ~ (hotlf "o)ci)) "o* = 

h 2 = (u H 2 + u 1 H 1 + u 2 H -h u 2 + (uottf HOtu + (ujfH ) (1) - - (7iJ%) m + 

+("otlf Ho)(2) - dottf "0)(2)) "o* 
= u H 2 u * = Vid C 4 

/i! vanishes as expected and h 2 simplifies to V, because u x , u 2 an< 3 H x vanish identically, and the product 
of two momentum-dependent functions contains no first- and second-order terms in l /c (equation ( 13.4D ). 
So far, we did not need to calculate one line explicitly to arrive at this result! The first three terms of the 
effective hamiltonian are obtained by sandwiching h to h 2 with 7i ref . 

Kffo = n re( h Q 7i ref = £id C 2 = \/m 2 + £ 2 id c2 

freffl = ^ref^l ^ref = 

Kin = n ie( h 2 7i ref = Vid C 2 

Finally, for h eS3 , we need to make some explicit computations and the first magnetic correction (third 
order in V0- There are three groups of terms which survive: 

heffS = ^iefh ^ref = ^ref ("3 #0 ~h U 3 + {u Q f c H 2 ) m ~ (M? U )(l) + ("ott?^o)(3) ~ (M? U o)(3)J "o* ^ref 
= : K«30 + fr e ff31 + Kff33 

The first two vanish when we project with 7i ref from left and right, because h eSQ = Eid c i is a scalar 
symbol, 

h ef{ 30 = ^ref (»3 H ~ h "3) u 0* ^ref = ^ref u 3 u 0* u H u 0* ^ref ~ ^ref K u 3 u 0* ^ref = °- 
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The second and third group of terms need to be calculated explicitly; since the details are arithmetically 
intricate, we have moved them to AppendixjE] (u j)^H 2 )n) — (h 2 ^ u o)(i) gives a gradient coupling to the 
potential, 

Kfm = Ttref (("oli? H 2 ) m - {h 2 f c u ) (1) ) n ref = -in iei {u , V} 7I ref 

= 9FfF 1 , r (V x VAg)-g. 
2£(E + mj 

The last term, h eS33 , contains the spin-orbit coupling: 

keff33 = TTref (Kf^o)(3) ~ ( h of "o)(3)) "o* ^ref = J B ljM n ref (%"o %j#0 ~ % ^0 %j"o) u 0* ^ref 
1 

= B ■ a 

2E 

Altogether, the effective dynamics up to errors of fourth order in 1 /c are given by 

ft eff = Eid c2 + iv id c2 + ~ ( — J— -(V X V A • a - -"-B • a") + tf(i/c<). (3.16) 
c A c" > \2E{E + m) 2E J 

The third-order correction is responsible for the spin dynamics and leads to the so-called T-BMT equa- 
tion. This result has been previously derived by Cordes IICor83ll and Teufel IITeu03l Section 4.1] under 
different hypothesis. We reiterate that the physical mechanism underlying the adiabatic decoupling in 
Teufel's work is different from the mechanism here. 

If we want to make this result rigorous, we will have to explicitly show that the construction of space- 
adiabatic perturbation theory still works when one replaces usual Weyl calculus with magnetic Weyl 
calculus. This has been the motivation for making the two-parameter expansion rigorous in the first 
place, but deserves a publication in its own right IIFL10II . 



A Equivalence of Weyl systems in both scalings 

Lemma A.l The adiabatic scaling and the usual scaling are related by the unitary U e , (U e <p)(x) := 
£- d/2 ip(^), ip G L 2 (R d ), i. e. we have 

Q = U e Q e U; 1 
P A , = [LIT 4 , U~ l . 

Proof Let y e L 2 (M d ). Then we have for Q e 

(17, Q e U^U^Xx) = (U E Q^){x) = e-^ (Q e <p) (f ) 
= s-^s^( x -)=Q(U e ip)(x). 

Similarly, we get for the momentum operators 

<x u^u^Xx) = (u e u A e> ^)M = «r* (ir^) (£) 

= e-^ (-i(V^)(f ) - ) V (f )) = (-isV x - M(Q)) (U eV )M. 

Hence the two scalings are unitarily equivalent. n 

Corollary A.2 The Weyl systems associated to the two scalings given by equations dl.lt and (II .2D are 
unitarily equivalent. 

Theorem A.3 The asymptotic two-parameter expansions of the magnetic Weyl products with respect to 
either scaling are given by the same terms order-by-order in e and A. 
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Proof To show that the asymptotic expansion of the product is the same, we have to revisit Theorem l2.10l 
(proof of equivalence of the two non-asymptotic product formulas product formulas) and translate the 
relevant formulas to the usual scaling. It suffices to show that the twister in both cases is the same 
function and thus the expansion has to be identical, too. We denote magnetic Weyl quantization with 
respect to the Weyl system in usual scaling, W A {X) := e Icr ( x >W' p e,;.M ; with Op^. For convenience of the 
reader, we will follow the notation in the proof of Theorem ll.ll as closely as possible. 
With a simple scaling argument, we get the composition rule for the Weyl system W A {X): 

WfrY) W^Z) = W^'Xey, r,) W^'\ez, Q 

= gi^sy.nUa.CSi nV«*(Q J Q+ e y,Q + ey + ez) W l/ '\ey + ez,r) + Q 

= e^^fi* A (Q, Q + ey,Q + ey + ez)W*(Y + Z) 

In Step 1 of the proof, we conclude from the composition law of the Weyl system (reformulated in the 
usual scaling), 

f dY {&; X f){Y) (K'S^Z - Y) e l 2 a[Y ' Z1 - 



Op^Cf)O^(g): 
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2d 



dZ 



that we need to find the operator kernel for 

n B sX (Q,Q + ey,Q + ez)W A (Z). 
If we apply this operator to a function e L 2 (K d ), we obtain 



(f2f A (Q, Q + ey,Q + ez)W A {Z) V ){v) = fi* A (v, v + ey,v + ez) e 



-i(v+«/2zK „-iVer A ([v,v+ez]) 



(p(v + ez) 



du e -Ku-/»>C e -i»-Ur\[u-ez,u]) n B^ u _ ez> u + ey> u ) 5 ( u _ ( v + ez)) 



duK{y,Z;u,v)^(u). 

To find the symbol associated to this object, we employ the Wigner transform adapted to observables in 
the usual scaling defined by 

W\y,^){X) :=e d {& a (tp,WftW))(-X) 

dye- iy ^e- iX, ' rA(lx - c/2y ' x+l!M) if*{x - e -y)i}>{x + |y). 
If we apply this to the integral kernel above, we get by the essentially the same calculation as before, 



{W A K{y,Z;;-)){X) = e a 



du e -'"-« e - lA/f r*C[*-'/™,*+'/™]) K(y, Z ;x- e -u,x + f u) 



e d e~ d 



,ia(X,Z) 



£ ue -iu-Z e -i></ e r A ([x->'/2u,x+e/2u]~) e -i(x-'/3u-*/2i)i) £ -i V«r*([3c-»/2 u-ez.x-'/i u]) . 

• fi* A (x - | u - ez, x - | u + ey, x - | u) 5(z + u) 

If we plug this into the remainder of the proof, we see that the twister term (after replacing z with y + z 
just as in Step 3) obtained here is identical to the one obtained in the adiabatic scaling, 

Hence the two expansions need to agree. □ 
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B Expansion of the twister 

Lemma B.l Assume B satisfies Assumption ^. 1\ Then we can expand j B around x to arbitrary order N in 
powers of e: 



Y 



((1 - (-l)" +1 )c - (1 - (-l) c )(rc + l)) y ; - ■•■y J -_ 1 z J - +i? N [r?](x, y,*) 

(B.l) 



n=1 |ai+|/3|=n-l 
N 

=:-Y i ^^n+RN[r B E Kx,y,z) 

71=1 



(B.2) 



Zn particular, the flux is of order e and the nth-order term is a sum of monomials in position of degree n+1 
andeachofthetermsisam^^{M d x ,^^{W i y xM d z ))function. The remainder is a ^ < ^ 00 (R d , ^ p ~(M d xR d )) 
function that is 0(e N+1 ) and can be explicitly written as a bounded function of x, y and z as well as N + 2 
factors of y and z. 



Proof We choose the transversal gauge to represent B, i. e. 



r 1 



A^x + a)- 



dsB; ; (x +sa)saj 



(B.3) 



and rewrite the flux integral into three line integrals over the edges of the triangle 
1 i- 1 



T e {x,y,z) = - 



dt [e (y, +z l )A l (x + e(t - i/ 2 )(y + z)) + 
-ey,A,(x + e(t - i/ 2 )y - §*) - ez { A, (x + |y + e(r - 1/2)2)] 
dss [-By(x + est(y +z)) (y, + z,) t(y } +z } )+ 



dt 

h 

+B l} (x + es{ty - §)) y t (ty ; - |) +B ;j (x + es(f + tz)) z, (f + tz ; )] 

All these terms have a prefactor of e which stems from the explicit expression of transversal gauge. We 
will now Taylor expand each of the three terms up to JV — 1th order around x (so that it is of Nth order 
in e). 

■ +1/2 r 1 



dt 



1/2 



dsB (j (x + est(y + z)) esCyy +z ; ) ; 



r i p+J/2 w-i e „+i 



ds 



df I] —s n+1 s~ n d Xh ■ ■ ■ d Xj B l]n+i {x) t" +1 Y\(y jm + Zj J + R im {x,y,z) 

V2 n=0 U ' m=l 

The remainder R 1N ; is of order JV + 1 in e, bounded in x and polynomially bounded in y and z. It is a 
sum of monomials in y and z of degree JV + 1. 



■(1 - zf- 1 3 T %-(x + £Tst(y + z)) est(y + z) 





f 1 


f 1 


p+1/2 


K iN>U,y, z ) = 


dT 


ds 


dt 










-1/2 



(JV-1) 



r.N + 1 



^ a! 



ri 



|a|=JV 



ds 



r 1 



dtst w+1 (y + z) a (y ; . + Zj .) 



dT (1 - Tf" 1 3«B (J (X + £TSt(y +z)) 
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The nth order term in e (the n — 1th term of the Taylor expansion) reads 



(n-l)l 



r i r +V2 
dss 

o J -1/2 



dt t» 3 X . • • • 3 Xjn _ By n (x) f](y Jm + Zj J 



ie n ny H i + (-i) n 



2 n! V 2 



— — -5 X ••■o> B u (x) V 



m 



The other factors can be calculated in the same fashion: 



(n-l)l 



r +V2 r i 
dt 

1/2 

sVl' " +2 
ni I 2 



o 



ds5»i-t"- 1 3 • • • d x . n _B ljn (x) Y\(ty jm - \z jm ) 



m=l 



hfn \ (-l) n - m + (-!)" 



m=0 



m+1 



The remainder is also of the correct order in e, contains JV+2 qs and a SS c £ co (R d , ¥™ 1 (E 1, xR d )) function 
as prefactor: 



R 2Nl ix,y,z) = s N+1 * ds dts (ty - |) a (ty, - f )• 

|a|=W Jo J -1/2 



\a\=N 

r 1 

dT (1 - zf- 1 (x + e T (sty - |)) 



The last term satisfies the same properties as 



+ 1/2 

dt 

-1/2 J 



r i 



d^T r- 

o (n-1)! 



m=l 



n! V 2 



^ f ! + (-!)"- 



mj n + 1 - m 



K 3N ; satisfies the same properties as R 1N ; and R 2 n!> 



B 3n; (w) = * N+1 J] ds dts (I + tz T (T + tz ;) 



a|=N 



a! 



ds 
o J 



f 1 



dT (1 - Tf" 1 3;S (; (x + £Ts(| + t Z )) . 



Put together, we obtain for the nth order term: 



1 e n n 



2 n! V 2 



m=0 v J 



l + (-l) n (-l)"- m + (-l) n l+(-l) n " m 

(ji # z ' 



n+l 

n+l 



m+1 



n + 1 — m 



^•••30,„ z ; m+1 --- z ;„ 



n! V. 2 J (n + l) 2 



^ 3;3/B„(x)y ;Zfc - 



|o|+|j9|=n-l 



(|a| +\] ((1 ~ ^-^'"'^X" + 1) - (1 - (-D n+1 X|a| + 1)) J V 
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The total remainder of the expansion reads 



R N [r B E ] = Ji (Rim -Rini) +z, (Rim-Rsm) e ^ M (K d ,<*?~(K d x M d )). 

In total, the remainder is a sum of monomials with bounded coefficients of degree JV + 2 while it is of 

□ 

B 



C Properties of derivatives of y 

For convenience, we give two theorems found in [IMP07] on the magnetic flux and its expontential 
which are needed to make the expansion rigorous: 

Lemma C.l If the magnetic field B (j , 1 < l,j < n, satisfies the usual conditions, then 

3 Xj J B = D jk (x, y, z) y k + E jk (x, y, z)z k 

d yj B E = D 'jk( x ' y> z ) yk + E ' jk ( x > y> z H 

9 *,rf = D'f k ix,y, Z )y k + E'; k {x,y,z)z k 
where the coefficients D ]k , E" k e S& <g°°(m. d x R d x R d ), 1 < j, k < d. 

Proof The corners of the flux triangles of F B found in ' IMP07' differ from those of y B e , but the proof 
carries over with trivial modifications. □ 

A direct consequence of this is the following simple corollary: 

Corollary C.2 If the magnetic field satisfies the usual conditions, then 

l^^^ 6- ^ 1 ^" 3 ! ^ CaftcC<^> + <^>) lal+ ' 6|+|C| ^ C a6c <^) |a|+|b|+|c| <^> |a|+|b|+|c| Va,b,ceN d , 
i. e. derivatives of e~ l ^ r "^ x ' y ^ are functions in y and z. 

D Existence of oscillatory integrals 

To derive the adiabatic expansion, we have to ensure the existence of two types of oscillatory integrals, 
one is relevant for the {n, k) term of the two-parameter expansion, the other is necessary to show exis- 
tence of remainders and the fcth term of the A expansion. 

Lemma D.l Let f e 5" p "\ p e [0, 1]. Then for all multiindices a, a e N d 



G(X) 



(27l) d 



&Ye i < x - r) y a r ) a {&- 1 f){Y)=({-idz) a {+id x Yf){X) (D.l) 



exists as an oscillatory integral and is in symbol class SP™ ' a ' p . 

Proof Since / is a function of tempered growth, we can consider it as an element of 5"'(M 2d ). Then, 
we can rewrite G as G = & a x a E, a & a where x and f are the multiplication operators initially defined on 
y(R 2d ) which are extended to tempered distributions by duality. Then for any t^> e 5^(IR 2d ), we have 

(G )¥ >) = (J^f a J^/,¥>) = = (/,(+i3 5 ) a (-iaj» 

= ((-id^n+id x ) a f,ip) 

where (•, •) denotes the usual duality bracket. 

Thus, the integral exists as an oscillatory integral. G = (— id^) a [+id x ) a f is also in the correct symbol 
class, namely zym-mp > and the lemma has been proven. □ 
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The next corollary is an immediate consequence and contains the relevant result for the term-by-term 
expansion of the magnetic product. 

Corollary D.2 Let f G 5" p "\ g e 5" p m2 , p e [0, 1] and a, a, b,/3 e Nq be arbitrary multiindices. Then for 
all functions B e ^Sf 00 (R d ) the oscillatory integral 



G(X) :-- 



dY 



exists, is in symbol class y™i+™2-W+\b\)p an dyields 



dZe , CT (x,y+z) B(x)j a T} a (jr -i /Xy)z b^ ( ^-i gXz) (D2) 



B(x) {{-id^\+id x Yf)(X) ({-id^\+id x f g)(X). 



(D.3) 



In the proof of Corollary ID.2I we have used that we could write the integrals as a product of two in- 
dependent integrals. There is, however, a second relevant type of oscillatory integral that cannot be 
'untangled.' Fortunately, we only need to ensure their existence and not evaluate them explicitly. Again, 
we will start with a simpler integral over only one phase space variable and then extend the ideas to the 
full integral in a corollary. 

Lemma D.3 Assume / e 5^% g e Sf™ 2 , p G [0,1], e e (0,1] and e [0,1]. Furthermore, let 
G z < e SB^ 00 ^*, ^(M. d y x )) be sucb thatforall c,c',c" e N d 

|3;3;'3/'G T ,(x,y, Z )| < C ccV ,«y) + (z ))' c l + l c 'l + l c "' 

holds for some finite constant C cc < c « > 0. Then for all a, a, b, [3 e Nq and t, t' e [0, 1] 
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2d 



dY 



dz gicr(X,y+Z) e iT§cr(y,Z) , 



G^y^yVzTm/XnC^gXZ) (D.4) 



exists as an oscillatory integral in $t>™-i +m 2 (|a|+|f>l)p_ j ne ma p ( T; T ') _> j^, j 5 continuous. 
Proof Let us rewrite the integral first, the result will serve as a definition for the oscillatory integral J TT / : 
1 
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4d 



dY dY dZ 



dZ ((+i3 jf ) a (-ia p ) a e i<70r - y - y) ) ((+i3 ? ) b (-i5 5 )' 3 e I<7(;f - z - z) )- 



•e-f CT ^G T ,(x,y, Z )/(Y)g(Z) 



(2ti) 



4d 









dY 


dY 


dZ 









dZ e" 7 ^" 7 - 7 ) e io-Cx-z,z) e iT|a(r,z). 
• G T ,(x,y, Z ) ((-i^) a (+i%) a /)(Y) ((-i3 e ) b (+ia 2 /g)(Z) 

By assumption d"d£f e 5^ m i-l a lp a s well as d^d^g e 5^2-1% anc i we see that it suffices to consider 
the case a = b = a = /3=0. In this particular case, we estimate all seminorms: let n,v eNl Then we 
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have to bound 



4d 



a+b+c=n 
a+(S=v 



a+b+c=n 
a+(3=v 



(271) 



(27l) 4d 



d7 



dY 



dZ 



d c x GAx,y,z)f{t)g{z) 



dY 



AY 



dZ 



d 2 e icr(x-y,y) e i(T(x-z,z) e ixf <7(y,z). 



3; G T / (x, y, z) 35 25f (?) 3/ 3^g(Z) 



2 



dy 



a+b+c=n J 
a+p=v 



dz 



from above by an integrable function. To do that, we insert powers of (y)~ 2 , (z)~ 2 , (tj)~ 2 and (Q~ 2 via 
the usual trick, e. g. (y)~ 2 (l — A IJ )e lr, y = e IT) ' y . To simplify notation, we set L y := 1 — A y ; L z , L v and 
are defined analogously. Then, we have for any N 1 ,N 2 ,K 1 ,K 2 e N 



a+b+c=n *- Z7T J 
a+/3=v 



dy 



drj 



dz 



• ^G T <(x,y,z) 8«d«f fx - f z, I - t,) 3 x b 3/g(x + f y, ? - £) 



Cn'ft' 



a+b+c=n 
a+/3=v 
|a'|<2N 1 ,|/3'|<2W 2 



dy 



dr; 



dz 



(y)- 2 ^(z)- 2iV2 3;G T ,(x,y,z)- 



■ d«d^ +a f (x - h z , ? _ r,) a» s p+p g ( x + ^ y> ? _ ? ) 



dy 



drj 



dz 



dfe'^e'f"* (y)- 2Nl (z)" 2N2 - 



a pa be a be 

a+b+c=n 
a+p=v 
\a'\<2N 1 ,\p'\<2N 2 
|o'|+|6'|+|c'|<2K a 
\a"\+\b"\+\c"\<2K 2 

■ (vr^ (Zr 2K2 ^(y)^b«{zW x d*'d?G,,{x,y >Z y 

■ d x a+a "d^ a 'f (x - f %x - r,) d x b+b 'd^ +p, g (x + f y, ? - c) (d.5) 

Here, the bounded functions t^ Na are defined by d x {x)~ 2N =: (x) _2N y> Na (x) for all JV e N , a e Nq, and 
the constants appearing in the sum are defined implicitly. We now estimate the absolute value of each of 
the terms in the integral in order to find N 1 ,N 2 ,K 1 ,K 2 e N large enough so that the right-hand side of 
the above consists of a finite sum of integrable functions. Using the assumptions on G T - and the standard 
estimate (£ — ri) m < (£) m {r))' m ', we can bound the integrand of the right-hand side of (1D.5D in absolute 
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2K l( y\-2K 2 . 



{r])- 2K HO 



value by 

y fabcap- , ,_2N 1 +|c|+|c'|+|c"|/_\-2iV 2 +|c|+|c'|+|c"| 

a+b+c=n 
a+p=v 
\a'\<2N l ,\p'\<2N 2 
|a'|+|b'|+|c'l<2Ki 
|a"|+|b"|+|c"|<2K 2 

. _ ^mi-CM + la'Dp ^ _ Qm 2 -(\p\ + \p" \)p 

< C (£) m i +m 2-l 1 'IP ^y^-2JVi+|n|+2K 1 +2K' 2 ^ z j-2JV 2 +|n|+2K 1 +2JC 2 + ^^-2K 2 + |m 2 -|/3|p| 

a+/3=r 

< C(£) mi+m2 " |rlp (y)- 2N i + l"l+ 2if i+ 2;f 2 ^ z ^-2N 2 + |n|+2K 1 +2if 2 ^j-^+KI+lvlp ^-2it 2 + |m 2 |+|r|p_ 

Choosing and K 2 such that —2Kj + \ mA + \ v\p < — d, j = 1,2, ensures integrability in 17 and £. Now 
that and K 2 are fixed, we choose JVj and JV 2 such that -2JV ; + \n\+2K 1 +2K 2 <—d and the right-hand 
side of the above is an integrable function in y, 17, z and ( which dominates the absolute value of 0D.5D . 
Thus, we have shown 



m 1 +m 2 — |v |p 



for all n,v e Nq and hence J TT ' exists in y™i +m 2 if the exponents of y, 17, z and £ in equation QD.4H all 
vanish, a = a= b = f3=0. Similarly, for general a, a, b,/3 e N^, we conclude 7 TT / e 5 ^n 1 +m 2 -(|a|+|6Dp_ 
Since the above bounds are uniform in z and z' , the continuity of (t, z') — > J TT / in the Frechet topology 
of 5^ m i+ m 2-([a|+|f>Dp follows from dominated convergence. n 

E Details of calculations in example 

For convenience of the reader, we present the derivation of h eii31 and h e{{ 33 in more detail. The calcu- 
lation simplifies tremendously when one discards blockoff diagonal terms as soon as possible. To that 
effect, we use that a ; and anticommute, /3 2 = id C 4 and {E, ■ a)(£ ■ a) = £ 2 id C 4. Furthermore, we 
introduce the spin operators 



Pi 



<jj 
CT, 



7 = 1,2,3. 



Then, we plug in the definition of 7i eff31 and keep only terms that are purely blockdiagonal, 



'eff31 



-iS Xl V 



L rc( 



^2E{E + m) 

?i(2£ + m) 



— id C 4 (E + m) id C 4+ 

2V2E 5 HE + m) 1/2 



+ 



2v / 2£ 5/2 (£ + m) 3 / : 



^ref- 



^2E{E + m) 

Writing out E = \f m 2 + E, 2 and using (V X V ■a){E,-a) = (V X V* • £)id C 4 + i(V x V A ■ p, we get 

1 



heffsi = -i9 Xl Vn lei 
i(V x V-?) 



m£, ^(2£ + m) 

' Trf" ld c 4 _ TrfTEH ^2? ld c 4 + 

477J 4£ (£ + mj z 



2£(£ + m) 



L rcf 



4£ 3 (£ + m) 2 
2£(£ + m) 



(2mE 2 + 2E{m 2 + E 2 ) - 2£ 2 (£ + m)) 7i ref - 



2£(£ + m) 



^ref(V x VA^)-p7T ref 
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Similarly, we can compute h e ff 33 , 

heff33 = TTref (("ott? W o)(3) ~ C^ottf "o)(3)) "o* "ref = ^ref (%"o % H ~ hS 5 C,"o) "o* "ref 

mix ., ?i(2£ + m) 



i; "ref 



2v / 2£ 5/2 (£ + m) 1 ^ c 2 v / 2£ 5/2 (£ + mfl 
ml j l } {2E + m) 



■ id r 4 + ■ 



2 72E 5 / 2 (£ + m) 1 / 2 u 2 V^ 2 (£ + m) 3 / : 
Once we use B l j l x 2; ,■ = and plug u * back in, we get the claim, 

mil It&E + m) 



(I ■ a)P 



a j j3 



a,/3 a.j + 



yj2E{E + m) 



1 Q JL re{ . 



^eff33 — 2 B 'J ^ref 



iB, 



■a,+ 



2v / 2E 5/2 (E + m) 1 / 2 ; ' 2v / 2£ 5/2 (£ + mf 2 

1 - ?! 



^/2£(£ + m) 



V2E 3 HE + mfl 2 



u Ttref 



8E 3 (E + m) 
i2B y £, ? J (2£ + m) 



^ref (-?;(2£ + m) a ; - (£ ■ a) - £,(2£ + m){l • a)a ; - + 2E 2 (£ + m^aj /3 7i ref 



8£ 3 (£ + m) 



1 1 

2£ 2E 
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